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The electroweak structures of the FCNC-induced rare and radiative K decays are presented. While 
the former decays are sensitive to short-distance physics and offer exceptional probes for new 
physics, the latter are dominated by long-distance physics. Even so, they should not be set aside 
since they constitute essential ingredients to control the hadronic uncertainties occurring for the 
rare decays. This is illustrated by systematically reviewing the phenomenological strategies cur- 
rently in use to cleanly predict the rare K decay rates in the Standard Model. 
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1. Introduction 

Flavor changing neutral currents (FCNC) are among our best windows into physics beyond 
the Standard Model (SM). Indeed, they are absent at tree level and thus entirely generated through 
quantum loops. This gives them a particular sensitivity to electroweak scale physics. Thus, if New 
Physics (NP) occurs at a scale not much higher, one can expect significant deviations []IJ|. 

Rare and radiative K decays permit to probe the quark-level s —>■ d FCNC transitions, a fun- 
damental step for reconstructing the flavor sector of a NP model [§]. However, the typical scale of 
K physics is well within the QCD non-perturbative regime, and the difficult task of dealing with 
hadronic effects has to be faced. For this, the main tool is Chiral Perturbation Theory (ChPT), the 
effective theory based on the symmetries of QCD which describes the interactions between the light 
pseudoscalar mesons. The price to pay for not having solved QCD is a series of phenomenological 
low-energy constants which have to be fixed from experiment. 

Here, our purpose is to review the strategies which are or could be used to resolve the hadronic 
uncertainties occurring in the rare K + — > n + vv, Kl — > 7T°vv, Kl — > 7t°£ + £~ and Kl — > A t+ A t_ 
decays, in particular with the help of radiative K decays. This will be done in Sees. 3, 4, and 5. 
Before that, in the next section, we describe the electroweak structures driving all these decays. 

2. Short and long-distance effects in rare and radiative K decays 

The s — > d FCNC transitions depicted in Table [I] contribute to both types of decays, whose final 
states can actually be identical. What distinguishes them is their sensitivity to short-distance (SD) 
physics. Rare decays are mostly induced by the c and f -quark contributions, and possibly by NP. 
By contrast, radiative decays are fully dominated by the long-distance (LD) w-quark contribution, 
which makes them rather insensitive to NP. They are always driven by photon penguins (barring 
pure bremsstrahlung-type processes), hence their name. Note that photons can be real or virtual, 
and some radiative modes with Dalitz pairs are actually more rare than rare decays. 

Standard Model electroweak processes: Which of the up-type quarks contribute the most, as 
indicated in Table [T[ is found by combining the behavior of the FCNC loop as a function of the 
quark mass [||, with the scalings of the CKM elements. Specifically, the Z penguin (W boxes are 
understood) is always SD. Its dominant piece comes from a quadratic breaking of SU (2)l, leading 
to a m 2 uct jM\ / behavior for the loop function. The y penguin does not suppress the light-quark 
contributions, hence is SD only if CP- violating since then lm(V* s V u d) = 0. Finally, the 77 penguin 
always gives CP-conserving LD contributions; its CP-violating part being strongly suppressed by 
the additional heavy quark propagator. Even though of higher electroweak order, this process can 
be competitive thanks to the LD enhancement and large CKM coefficients. 

The decays in Table |l| are written in terms of the CP-eigenstates K{ and K2', the neutral mass 
eigenstates being Kl ~ K2 + sKi and K$ ~ K\ + £#2- Neutral K decays thus proceed either directly, 
for example Kl — > K2I— ► X], or indirectly, Kl — > £Ki{—> X]. The indirect contribution is always 
proportional to the small mixing parameter Sk = 2.228(11) x 10~ 3 , but is not always negligible 
since the direct contribution can be significantly suppressed. 

The K + decays are not included in Table [j], since CP is no longer at play. Instead, K + decays 
receive contributions from the mechanisms of both columns of Table [j]. For instance, K + —>■ 7l + vv 
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CP-violating CP-conserving 

t, c (+ small u) quarks: 
K\ -> 7T°VV 

Ki -» TZ^+r 

w quark dominates: 
JTi -► TZ^+r 

£ — > 7T^7 



m quark dominates: 

Ki,2 -> 77, ^1,2 -> 7T°77 
^i,2 -» 7r £+r 

*1,2 -> 

Table 1: Electroweak processes relevant for rare and radiative K decays, with the dominant contributions 
for each mode (7 in the final state stands for both a real photon or an Dalitz pair, and £ = e,)J,). 

has a "CP-violating" Z-penguin contribution dominated by the t and a "CP-conserving" Z-penguin 
contribution with t and c quarks. As can be seen in the table, for all other K + decays, there is a 
dominant CP-conserving w-quark 7 or 77 penguin, so all these decays are radiative. 

Windows for New Physics: In summary, the only modes giving us access to SD physics are: 
K + —> K + vv and K L — > 7T°VV (K s —> 7t°vv is difficult to measure), induced by the Z penguin, K L — > 
7tPe + e~ and Ki — > 7r°/i + ^t _ , which receive contributions from the three electroweak processes, and 
Kl — > ' f° r which the 7 penguin is absent (the helicity suppression makes ^ — > e + e~ very 

small). These modes are naturally the most sensitive to NP effects (see e.g. ||l|, ||, ||, ^] for the 
impacts of supersymmetry). Further, taken in combination, they offer a powerful discriminating 
tool thanks to their different sensitivities to the underlying electroweak processes. 

In principle, the SD parts can also be accessed through CP asymmetries. As will be briefly 
discussed in Sec. 4, this is especially clean for the K + modes, for which the direct CP-asymmetries 
typically arise through the interferences between the "CP-violating" and "CP-conserving" pieces 
of Table 1. For the neutral modes, CP asymmetries are often driven by the mixing and essentially 
proportional to Ek, with small corrections from direct CP-violating effects. 

3. Long-distance effects in K + — > tt + vv and K L — > 7t°vv 

Matrix elements: Even if the t and c-quark Z penguins are said to be SD, there are hadronic 
effects to be dealt with. Indeed, these SD contributions are encoded as effective four-fermion 
operators, of which we need the matrix elements between hadron states to get down to physical 
observables (see Fig. [ija). In the SM, there is only one dimension-six operator [|l|]: 

Qeff = sf (1 - 75)^® VYn (1 - 7s) v • (3. 1) 

To obtain (nvv\Q e ff\K), we need the matrix element of the vector current, (7i\sy^d\K), which is 
described in terms of the vector and scalar form-factors (only the former is needed for massless 




S W ± d 



?-quark dominates: 
K 2 -> 7T°VV 



u,c,tP" 



Only t & c quarks: 
K — > Tiny 




Only t & c quarks: 
negligible 

( ~ l/m c , t ) 
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Figure 1: LD effects in K — > nvv. a: The matrix elements of Q e ff and their relation to those of the charged 
current, b: The up-quark Z penguin and associated LD processes. 



neutrinos). These form-factors are directly related to those of the charged current decays K — > %l\ 
(so-called K^), induced by the Fermi interaction: 

QFermi = Sf (1 - ft) If ® Vty (1 - ft) I . (3.2) 

In the isospin limit, (7r|^ r y^«|A') = (7rj.s r y' , ci|A'), up to simple Clebsch-Gordan coefficients. 

Isospin-breaking effects have to be accounted for given the percent level of precision we are 
after. For that, three very clean theoretical ratios are used, involving the vector form-factors f'/(q 2 ) 
of the (n' \sj^u\Kj) and (n'\sY^d\K J ) matrix elements (q 2 is the momentum transfer), as well as 
their slopes X'/ defined as the coefficients of their Taylor expansions around q 2 = [| 



1.0000(2), r K 



/| + (Q) 



1.0015(7), r x 



-,00 



X 



h0,0H 



: 0.990(5) . (3.3) 



Using the results on the form-factors of the Flavianet fit to Kfj data the phase-space integrated 
matrix elements K v ~ / d^^\(jivv\Q e ff\K)\ 2 (see ^ for the precise definition) are obtained with 
an impressive precision: = 0.5168(25) x 10~ 10 and = 2.190(18) x 10~ 10 , with the errors 
coming mostly from data. Consequently, the errors from the LD effects on the t and c-quark 
contributions to K — > nvv amount to only a few percents of the total error [jl|]. 

Up-quark Z penguin: The M-quark contribution to the Z penguin is a priori suppressed compared 
to that of the c-quark by the tiny m 2 /m 2 ratio. However, since the w-quark is dynamical, this 
translates as a significant tff{m 2 K )/m 2 ~ 10% correction. At this level, the other LD contributions 
arising from the W box or to the AS = 1 ® AS = 1 tree-level processes also have to be included [ 10 , 
11], see Fig. [l]b. All these LD contributions are CP-conserving, hence do not occur for K L — > n°vv. 

The computation is done in ChPT by including the hadronized AS = 1 three-flavor effective 
operators Q\ to Q(, (in terms of the weak low-energy constants Gs and G27, fixed from K — > %% and 
thus accounting for the Al = 1/2 enhancement), together with the hadronized effective interactions 
ud,us —>W—>£v and iiu,dd,ss — > Z — > vv. A complication is that the standard prescription to 
include the neutral current within ChPT by promoting ordinary derivatives to covariant ones does 
not work in presence of AS = 1 interactions, because the GIM mechanism is missed. As a result, the 
contributions from (nvv\Q e ff\K) described in the previous section are entangled with the genuine 
w-quark Z penguin effects (in other words, ChPT regenerates the meson loop of Fig. [l^). 
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Figure 2: LD effects in — > n IJ £ + £ . a: The matrix elements of $ e ft and their relation to those of the 
charged current (DCPV). b and c: Meson representations for the 7 and 77 penguins (ICPV and CPC, resp.). 

A simple procedure to enforce the GIM mechanism is to require the absence of a local Kl — > Z 



coupling [JlOp , so that the meson loops have the expected w-quark Z penguin electroweak structure. 
Many counterterms occur, some of them related to K + — > 71+7* [— > Allowing for a conser- 

vative uncertainty for those which are unknown, the correction from the light degrees of freedom 
is 8P U = 0.04(2) to be compared to the c-quark contribution F c = 0.372(15) [JTJ, 

At this point, the other correction of order m\ / m c can be evaluated. It arises from dimension 
eight operators induced by the c-quark, i.e. four-fermion operators with two derivatives (each 
counting as G(mx) given the momentum scale involved). The problem reduces to the evaluation 
of the matrix elements of these operators. Though this cannot be done exactly, an approximate 
matching with the non-local w-quark Z penguin computed in ChPT is possible, making use of the 
CKM unitarity (the f-quark piece is negligible). It turns out that this correction 8P C is smaller 
than 8P U as it misses the AI = 1/2 enhancement. The combined correction is then simply 8P UC = 
0.04(2) and amount to a 6% increase of the K + — > K + vv rate [JTTJ| . 

4. Long-distance effects in K L — > n°e + e~ and K L — > 

These decays receive three contributions (see e.g. [||]). The first piece is the so-called direct 
CP-violation (DCPV), induced by the Z and 7 penguins with t and c quarks (see Table |I|). The 
second piece is the indirect CP-violation (ICPV), which proceeds through the small K\ component 
of the Ki. The — > n°l + £~ decay is CP-conserving and dominated by the LD w-quark 7 penguin. 
Both the DCPV and ICPV can produce the lepton pair in a 1 state, and thus interfere. Finally, 
the third piece conserves CP (CPC) , and proceeds through a purely LD 77 penguin. 

Each of these three pieces involves LD physics at some stage. Let us describe the strategies 
used to bring them under control. 

Matrix elements for DCPV: The situation is similar to that for K —> ttvv, with the SD Z and 7 
penguins with t and c quarks encoded into two effective dimension-six operators 

Q v eff =sf(l-y 5 )d® ly^i , Qfc f f = sf{\-y 5 )d® ly^i . 



(4.1) 



The hadronic matrix elements are the same as for K — ► Ttvv, though now the scalar form-factor 
contributes to (71°h + h~\Q^^\Kl) (its impact for 7t e + e~ is suppressed by the electron mass). 
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Isospin-breaking effects can be controlled as before (Fig. ^?), since the scalar /q (q 2 ) and vector 
f'-liq 1 ) form-factors are equal at q 2 = 0, while the scalar slope is not very precisely known but 



V,A 



has a small impact. Using the Flavianet fit the phase-space integrated matrix elements K { 
J 1 (7l Q £ + £~ I Q^ff \Kl) \ 2 are found with a few per mil precision (see [0] for details). This is more 
than enough given the larger errors on the other two contributions to Kl — > 7l Q £ + £~. 

ICPV contribution: The process K L -» sK { [-> 7i°£ + r] is related to K s « #1 -> 7r°^+£", which 
is dominated by LD physics (Fig. ||&). The ChPT analysis (supplemented with dispersion relations) 
of Ref . [ |l3| ] has shown that meson loops are subleading and the amplitude is essentially constant, 
proportional to a linear combination of counterterms denoted as- Using the experimental Ks — > 
K°e + e~ and Ks —> 7r°/^ + /i _ rates [14] gives \as\ = 1.2(2). The sign of as is inaccessible since 
loops are small and these rates are simply proportional to a 2 s . 

This 20% uncertainty on as is the main source of error for the Kl —> 7t°£ + £~ rates & Further, 
the sign of as determines that of the ICPV and DCPV interference, thus unknown at present. To 
improve our knowledge of a$, there are two possible approaches besides better measurements of 
Ks — > 7i°£ + £~. First, as can be related to the counterterm a + occurring for K + — > 7t + £ + £~, though 
in a somewhat model-dependent way [15, 1(3]. Second, the decay Kl — > n°n £ + £^ depends on the 



same counterterm as, and has a non-negligible loop contribution [|17[]. Though its branching ratio 
is small, of (ff(l0~ 9 ) for £ = e, its measurement would give us both the size and sign of as- 

Finally, it should be mentioned that, if statistics is sufficient, the rare decay Kl — > 7T ( V + A t 
itself could fix the sign of as through the observation of a forward-backward asymmetry [g], which 
is predicted to be of about sign(as) x 15%. 

CPC contribution: The final contribution is induced by the 77 penguin (Fig. ||c), and is purely 
LD. In ChPT, it is represented by ^ and A' ± loops, is finite at @{p A ), and produces the lepton pair 
in a ++ state only. It is helicity suppressed and thus relevant only for Kl — > 7T /x + /i _ . 

To get a handle on higher order effects, the experimental information on the Kl — > 71° 77 mode 
can be used (as well as on Ks — > 77, which involves the same ^(/? 6 ) counterterms). For the muon 
mode, the Kl — > 7T°77rate permits to partially account for &(p 6 ) effects, bringing the CPC piece 



under control to within 30% [18]. For the electron mode, the absence of the 77 contribution is 
confirmed by the photon energy spectrum in Kl — * % yy. Indeed, to escape the helicity suppression, 
the e + e~ pair is necessarily in a 2 ++ tensor state. This requires a significant production of the two 



photons also in a 2 ++ state, which would have had a clear signature on the spectrum 015Q. 
Direct CP-asymmetries: Radiative K + decays are LD-dominated, yet the asymmetries 

, , ,v r(K + ->x + )-r(K- ->x-) 

A CF (K^X+) = ^-^L-^— -J (4.2) 



are sensitive to SD physics. Possible channels are for example X + = K + £ + £ [ |13[ ] or X + = 
7T + 7r°7 [19]. However, these asymmetries are significantly suppressed in the SM by the smallness 



of the SD parts relative to the LD contributions, and should not exceed about 10~ 4 . Phase-space 
asymmetries, e.g. the Kl,s —> X M + M forward-backward asymmetry or the Kl — ► 7i + 7t~Y* angu- 
lar asymmetry, could also be interesting even though one needs to disentangle indirect and direct 
CP violation [§, 0]- 
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Figure 3: LD effects in Kl — > . a: The matrix element of 6^// an d its relation to that of the charged 

current, b: Meson representation for the 77 penguin, c: The pole amplitudes sensitive to Gg. 



5. Long-distance effects in K L — > 



At first sight, the absence of the 7 penguin appears to promise a relatively simple structure, 
and thereby a good control on the hadronic uncertainties. Indeed, the Z penguin is as always 
dominated by SD, and the decay being CP-conserving, it receives both t and c-quark contributions 
(see Table [j] and Fig. ^p). However, this picture is upset by the large contribution from the 77 
penguin (Fig. |3|fo), which gives rise to three problems to be addressed before being able to use the 
measurement B(Kl —> /i + ;U~) = 6.87(11) • 10~ 9 to constrain SD physics. 

First, the absorptive part of the 77 loop, precisely known from Kl — > 77, contributes for 6.7 • 
10~ 9 and nearly saturates the branching ratio. This clearly makes the current experimental precision 
on B(Kl — > pL + pL~) insufficient to access to the interesting SD physics. Second, the 77 loop is 
divergent in ChPT, in stark contrast to Kl — ► n°£ + £~ where it is finite at leading order. This means 
that unknown counterterms occur and the dispersive part of the 77 loop cannot be predicted within 
ChPT. Finally, this dispersive part interferes with the SD piece, but the sign is unknown. 

Though a definitive answer is still lacking, let us detail the current and possible future strategies 
designed to bring this mode under control. 

Matrix element: The effective operators are the same as for Kl — > n Q £ + £~, but with different 
SD coefficients since the 7 penguin is absent p0[]. Only the (n + n~ \Q^^\Kl) matrix element is 
non-zero, and is again related to that of the charged current (jJ. + v^ \QFermi\K + ) which, by PCAC, is 
given in terms of the K + decay constant (Fig. ^a). Isospin-breaking effects could be included but 
are not currently needed given the uncertainties on the other contributions. 

Sign of the 77 contribution: To get this sign, we need that of the Kl — > 77 amplitude. However, 
its evaluation is rather tricky because it vanishes at leading order in ChPT. This cancellation occurs 
because the operator driving Kl — > 77 is Q\ = sy^{\ — ys)d ® Sy^(l — Ys)u, and in a pole model 
(Fig. |||c), there is no way to get an uu state as a linear combination of 71° and T]s states [21]. The 



same mechanism is at play for AMls or for the non-local magnetic contribution to Kl — > K + K^y. 

To circumvent this problem and consistently account for the class of higher order corrections 
induced by the singlet Tji in ChPT, the best strategy is to first go to U(3) ChPT, where the pole 
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NP probe Rare decays Decays to be used to control LD effects 
K — > 71VV Kft (matrix elements) 

K + — > n + £ + £~ (up-quark Z penguin) 



X 



SZ,y,H k l -> 7i°£+r Kb (matrix elements, DCPV) 

^ ' K s ^n Q £ + r, K+ ^n+£ + £-, K L ^n°7i°£ + e- (ICPV) 

^ Kl ► 7r°77, -» 77 (CPC) 

j J (also accessible with ^+ — > 7T+£+i?~, ^+ — > n + n°y, ... CP-asymmetries) 



X 



Z H Kl —> n + n K12 (matrix elements) 

Kl — > 77, #s — ► 7T°77, #+ — > 71+77 (interference sign) 

-» 7^+r , #l -> £ + £~l' + t'~ (dispersive 77) 



Table 2: Sensitivities of the rare decays to NP entering Z, 7 and non-standard Higgs penguins [[[J |^, ||, 0]. 
The last column lists the modes used to control their hadronic uncertainties. 



amplitudes no longer vanish, and then project back onto standard SU (3) ChPT. Doing so for Kl — > 
77, the amplitude turns out to be proportional to a new weak low-energy constant G|, of which we 
now need the sign. To this end, one can either make some generic assumptions on the behavior of 
the <2i to <26 operators in the non-perturbative regime, or constrain G\ experimentally from those 
pole amplitudes which do not vanish at leading order in SU (3) ChPT, for example K$ — > ^ 77 or 
K + — > 71+77 (Fig- 11^)' see R e f- [0] f° r more details. 

Size of the 77 contribution: The problem reduces to that of computing the value of a low-energy 
constant. This is very difficult as one must deal with hadronic physics above the ChPT scale but 
still well within the non-perturbative regime. The standard approach is to parametrize the form- 
factor for the vertex ^z, — > 7*7*, and then use this form-factor to compute the 77 loop. For moderate 
photon virtuality, one relies on Vector Meson Dominance (VMD) or large Nc arguments (Fig. |]fr), 
while at high virtuality, an approximate matching with the partonic u and c-quark 77 penguin is 



performed (see e.g. [ ]22[ ] and references there). 

The VMD ansatz introduces a number of free parameters which have to be fixed from the Kl — > 
je + e~, Kl — ► 7jU+jU~ and Kl — > e + e~ £i+£i~ differential rates. This can only be done partially 
at present, and improving the experimental accuracy would be welcome, especially on the latter 
mode. Kinematically, Kl — > jU+ jU~ jU+ /i~ would be even better but its 10~ 13 branching ratio [23] 



presumably forbids any study of the differential rate. 

Altogether, the dispersive 77 contribution is similar in size to the SD part, but is poorly known 



and thus the main source of uncertainty on the Kl — > jU+|U rate in the SM [22] 



6. Conclusion 

The rare decays K — > 7TVV, Kl — ► 7t°£ + £~, and to some extent Kl — > jU+jU~ offer unique 
windows into 5 — > d transitions thanks to their sensitivity to SD physics, and thereby to possible 
NP effects. In addition, each mode can probe specific electroweak processes, and thus, when taken 
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in combination, these decays form a powerful discriminating tool (see Table Q). Even if NP is first 
discovered at the LHC, they will remain essential to uncover its flavor structures. 

In this program, radiative decays have an important supporting role. Being dominated by 
LD physics, they are mostly insensitive to NP. This makes them appropriate to reliably fix the 
hadronic quantities originating from residual LD contributions in rare K decays (see Table ||). The 
exceptional cleanness of the rare modes would not be possible without experimental information on 
the radiative ones (as well as on the Kb decays). Therefore, besides the golden rare decay modes, 
radiative decays should be an integral part of the experimental programs under development 
and will ensure a rich and fruitful K physics phenomenology for the years to come. 
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